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ABSTRACT 
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O: Aims. The physics of star formation is an important issue of Galactic evolution. Most stars are formed in high density 

environments (n > 10 4 cm -3 ) emitting lines of diverse molecular transitions. In the present part of our survey we search 
for ammonia emitters in the Aquila rift complex which trace the densest regions of molecular clouds. 
Methods. From a CO survey carried out with the Delingha 14-m telescope we selected ~ 150 targets for observations in 
other molecular lines. Here we describe the mapping observations in the NH3(1,1) and (2,2) inversion lines of the first 
S—j . 49 sources performed with the Effelsberg 100-m telescope. 

Results. The NHa(l,l) and (2,2) emission lines are detected in 12 and 7 sources, respectively. Among the newly discovered 
NH3 sources, our sample includes the following well-known clouds: the starless core L694-2, the Serpens cloud Cluster 
B, the Serpens dark cloud L572, the filamentary dark cloud L673, the isolated protostellar source B335, and the complex 
star-forming region Serpens South. Angular sizes between 40" and 80" (~ 0.04 — 0.08 pc) are observed for compact 
starless cores but as large as 9' (~ 0.5 pc) for filamentary dark clouds. The measured kinetic temperatures of the clouds 
lie between 9 K and 18 K. From NH3 excitation temperatures of 3 — 8 K we determine H2 densities with typical values 
of ~ (0.4 — 4) x 10 4 cm -3 . The masses of the mapped cores range between ~ 0.05 and ~ 0.5M Q . The relative ammonia 
abundance, X — [NHa]/[H2], varies from 1 x 10~ 7 to 5 x 10~ 7 with the mean (X) = (2.7±0.6) x 10~ 7 (estimated from 
spatially resolved cores assuming the filling factor 77 = 1). In two clouds, we observe kinematically split NH3 profiles 
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separated by ~ 1 km s _1 . The splitting is most likely due to bipolar molecular outflows for one of which we determine 
an acceleration of V < 0.03 km s _1 yr _1 . A starless core with significant rotational energy is found to have a higher 
kinetic temperature than the other ones which is probably caused by magnetic energy dissipation. 



X 



Key words. ISM: clouds — ISM: molecules — ISM: kinematics and dynamics — Radio lines: ISM — Line: profiles 
Techniques: spectroscopic 



1. Introduction tending up to latitudes of at least +10° and down to 

-20° (e.g., Dame & Thaddeus 1985; Kawamuraet al. 1999; 

One of the central questions of stellar evolution is about the Dobashi et al. 2005). 

processes defining the initial mass of a protostellar object. The distance to the Aquila rift was est i ma ted from opti- 
It is suggested that more than 70% of stars are formed in Cjd dat& &nd interstellar extinction, and from parallax mea- 
dusters embedded within giant molecular clouds (Lada k surements _ A sudden appea rance of reddened stars is ob- 
Lada 2003) and that the varying environment may affect seryed ftt D = 225 ± 5g pc &nd the doud distribution may 
the initial mass function (IMF). have & thickness of „ 80 pc ( Strai z ys e t al. 2003). This es- 
In the present survey we search for the star-formation timate is in line with the distance to the Serpens- Aquila 
regions in the Aquila giant molecular cloud. Starting region found from interstellar extinction by Knude (2011), 
from Cygnus, in the direction of the Galactic center, the D — 203 ± 7 pc, but significantly smaller than the trigono- 
Milky Way appears to be split into two branches. The metric parallax value 415 ± 25 pc deduced by Dzib et al. 
southern branch runs through Cygnus, Vulpecula, Sagitta, (2010) from the Very Long Baseline Array observations of 
Aquila, and Scutum, entering the Galactic central bulge the young stellar object EC 95 in the Serpens cloud. Since 
in Sagittarius. The northern branch crosses Vulpecula and the discrepancy in the distance estimations may be caused 
Aquila and disappears in the northern part of the Serpens by a projection effect, we will use hereafter the weighted 
Cauda and Ophiuchus constellations. The separating dark harmonic mean of the first two measurements, both based 
band of obscuring dust is commonly called the Aquila on large-size data sets, D = 203 ± 12 pc, if not otherwise 
rift, but absorbing clouds are found not only next to the specified. Thus, comparing to other close molecular corn- 
Galactic plane: many dark clouds are located in areas ex- plexes, the distance to the Aquila rift is similar to that 
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of the Perseus clouds (<~ 250 pc; Rosolowsky et al. 2008) 
and lies between those of the Taurus and the Pipe Nebula 
(~ 140 pc; Rathborne et al. 2008) and the Auriga (~ 300 
pc) molecular clouds (Heiderman et al. 2010; Cieza et al. 
2012). 

Recently the Herschel Gould Belt Survey reveals in the 
Aquila rift cloud complex a large number ( > 500) of star- 
less cores (e.g., Bontemps et al. 2010; Konyves et al. 2010; 
Mcn'shchikov et al. 2010). Starting our investigations of the 
Aquila clouds we have carried out in 2010 an extended sur- 
vey of CO and extinction peaks using the Delingha 14-m 
telescope (Zuo et al. 2004). The targets were selected from 
the CO data of Kawamura et al. (1999, 2001) and from the 
optical data of Dobashi et al. (2005), limited in the latter 
case by the areas with A v > 6 m . We observed and de- 
tected 23 Kawamura and 47 Dobashi clouds at A <~ 3 mm. 
However, CO and 13 CO trace not only the dense gas but 
also outflows and low density regions; in addition, CO may 
be saturation broadened. It is well known that the carbon- 
chain molecules disappear from the gas-phase in the cen- 
tral regions because of freeze-out onto dust grains (e.g., 
Tafalla et al. 2004), i.e., C-bearing molecules are usually 
distributed in the outer parts of the cores. That is why CO 
and 13 CO are not good indicators for the dense molecular 
gas. Contrary to this, N-bearing molecules such as ammo- 
nia concentrate in the inner cores, where the gas density 
approaches 10 5 cm -3 . Ammonia is still observable in the 
gas-phase since it is resistant to depletion onto dust grains 
(e.g., Bergin & Langer 1997). 

Based on the CO survey we composed a list of <~ 150 
targets where not only CO but also 13 CO is quite strong, 
thus indicating large column densities and hence promising 
conditions to detect dense cores. The objective of the cur- 
rent work performed with the Effelsberg 100-m telescope is 
to observe these targets in other molecular lines, in the first 
place in NH3, in order to confirm the target identifications 
as dense cores and to determine their physical properties. 

In the present paper we describe the first 49 sources 
from our survey Observations are described in Sect. 2. 
Individual molecular cores are considered in Sect. 3. The 
results obtained are summarized in Sect. 4, and computa- 
tional methods are outlined in Appendix A. Appendix B 
contains NBi3(l,l) and (2,2) spectra detected toward the 
ammonia peaks, tables with derived physical parameters, 
and maps showing distributions of these parameters within 
the three most prominent cores. 



2. Observations 

The ammonia observations were carried out with the 
Effelsberg 100-m telescope 1 in 4 observing sessions between 
March 26 and 30, 2011, and in 2 sessions on January 12 
and 13, 2012. The NH3 lines were measured with a K- 
band HEMT (high electron mobility transistor) dual chan- 
nel receiver, yielding spectra with full width to half power 
(FWHP) resolution of 40" in two orthogonally oriented lin- 
ear polarizations at the rest frequency of the (J, K) = (1,1) 
and (2,2) lines (23694.496 MHz and 23722.631 MHz, respec- 
tively). Averaging the emission from both channels, the typ- 



1 The 100-m telescope at Effelsberg/Germany is operated by 
the Max-Planck-Institut fur Radioastronomie on behalf of the 
Max-Planck-Gesellschaft (MPG). 



ical system temperature (receiver noise and atmosphere) is 
100 K on a main beam brightness temperature scale. 

In 2011, the measurements were obtained in frequency 
switching mode using a frequency throw of ±2.5 MHz. 
The backend was an FFTS (fast Fourier transform spec- 
trometer), operated with a bandwidth of 100 MHz, provid- 
ing simultaneously 16 384 channels for each polarization. 
For some offsets we also used the minimum bandwidth of 
20 MHz for the NH 3 (1,1) line. The resulting channel widths 
were 0.077 km s _1 and 0.015 km s _1 , respectively. We note, 
however, that the true velocity resolution was about 1.6 
times as large. 

In 2012, the measurements were obtained with the back- 
end XFFTS (extended bandwidth FFTS), operated with 
100 MHz and 2 GHz bandwidths, providing 32 768 channels 
for each polarization. The resulting channel widths were 
0.038 km s _1 and 0.772 km s _1 , respectively, but the true 
velocity resolution is 1.16 times larger (Klein et al. 2012). 

The pointing was checked every hour by continuum 
cross scans of nearby continuum sources. The pointing ac- 
curacy was better than 5". The spectral line data were 
calibrated by means of continuum sources with known 
flux density. We mainly used NGC7027 (Ott et al. 1994), 
3C123 (Baars et al. 1977), and G29.96-0.02 (Churchwell 
et al. 1990). These calibration sources were used to estab- 
lish a main beam brightness temperature scale, T MB . Since 
the main beam size (40") is smaller than most core radii 
(> 50") of our Aquila objects, the ammonia emission cou- 
ples well to the main beam and, thus, the T MB scale is appro- 
priate. Compensations for differences in elevation between 
the calibrator and the source were < 20% and have not 
been taken into account. Similar uncertainties of the main 
beam brightness temperature were found from comparison 
of spectra toward the same position taken at different dates. 

3. Results 

The detected NH3 spectra were treated in the way de- 
scribed in Appendix A. We note that the radial velocity, 
Ksr, the linewidth Av (FWHM), the optical depths, m 
and t 2 2, the integrated ammonia emission, jT MB dv, and 
the kinetic temperature, Tki n , are well determined phys- 
ical parameters, whereas the excitation temperature, T ox , 
the ammonia column density, A r (NH 3 ), and the gas density, 
n(H2) are less certain. The optical depth and kinetic tem- 
perature are parameters independent of calibration errors, 
since T^ n depends only on the intensity ratio of the simul- 
taneously observed NH3(1,1) and (2,2) lines, and t depends 
on the relative intensities of the various hyperfine (hf) com- 
ponents. On the other hand, the excitation temperature T ex 
depends on the ammonia line brightness and is sensitive to 
calibration errors (Lemme et al. 1996). Besides, Eq.(A.lO) 
shows that T ex depends on the filling factor 77 which is not 
known for unresolved cores. 

A typical error of the optical depth determined formally 
from the covariance matrix at the minimum of \ 2 is 15%- 
20% (Appendix A). The error of the integrated ammonia 
emission depends on the noise level and for an rms ~ 0.2 K 
per 0.077 km s^ 1 channel it is <~ 0.3 K km s _1 . The preci- 
sion with which the Gaussian line center can be estimated 
is given by (e.g., Landman et al. 1982) 

i Tins 

g " (2 7 rln2)V4 ^VA^, W 
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where A c h and Av are the channel and the line widths, 
respectively. At rms ~0.2K and T MB ^ 2 K (see below) it 
gives ex y <~ 0.005 km s _1 for our high resolution settings. 
The same order of magnitude errors (i.e., from a few to ten 
m s _1 ) of the line centers and linewidths are obtained from 
the covariance matrix. 

The results of our analysis are presented in Tables 1-2 
and B.1-B.4. Table 1 lists the regions where NH3 line emis- 
sion was detected either in both (1,1) and (2,2) transitions 
or only in (1,1). In case of non-detections (Table 2), the rms 
T MB values per channel width are given. In total, NH 3 lines 
are detected in 12 sources out of 49 dense cores. In Table 1, 
we list the radial velocities, Vl SR , the observed linewidths 
Av (FWHM), and the main-beam-brightness temperature, 
T MB , at the maximum of the main group of hf components. 
In Col. 2, different peaks resolved in the sources are labeled 
by Greek letters (a for the most intense NH3 emission) 
which are also indicated on color maps (Figs. 1-3, 5, and 
7-10). The physical parameters calculated for each individ- 
ual source are given in Table 1 and in the Appendix B, 
Tables B.1-B.4. The results obtained for the three most 
prominent sources Do279P6, Do279P12, and SS3 are shown 
in addition as maps (Figs. B.10-B.12). Below we describe 
the individual sources identified as dense molecular cores. 
Statistical analysis of the full sample of Aquila cores will 
be reported after the survey will be completed in 2013. 

3.1. Kawamura 01 (KaOl) 

The starless cloud Kawamura 01 (KaOl herein) belongs to 
the dark cloud L694 classified as opacity class 6 by Lynds 
(1962). This is an isolated dense core located in the sky 
close (<~ 4°) to the protostar B335 (Sect. 3.6). No embed- 
ded IRAS luminosity source has been found (Harvey et al. 
2002; Harvey et al. 2003a). Lee et al. (1999) found three 
dark spots L694-1, L694-2, and L694-3 in the Digital Sky 
Survey image of L694. The second spot is associated with 
KaOl. Using Wolf diagrams, Kawamura et al. (2001) deter- 
mined the distance to L694-2 as D m 230 pc. 

We mapped KaOl for the first time in the NH 3 (1,1) and 
(2,2) lines with a spacing of 40" at 29 positions marked by 
crosses in Fig. la. It is seen that the NH3 emission arises 
from a very compact region. The fourth contour level in 
Fig. la corresponds to the half-peak of the integrated emis- 
sion, fT MB dv. This allows us to calculate the apparent geo- 
metrical mean diameter (FWHP) d = Vab, where a w 115" 
and b rts 65" (deconvolved) are the major and minor axis of 
this annulus, respectively, and a beam ammonia core size of 
approximately 0.08 pc — the typical value for dense molec- 
ular cores in the Taurus molecular cloud (e.g., Benson & 
Myers 1989). 

At the core center, we obtained two spectra with low 
(FWHM = 0.123 km s" 1 ) and high (FWHM = 0.024 
km s _1 ) spectral resolutions. The high resolution spectra 
of the NH 3 (1,1) and (2,2) lines are shown in Fig. B.l along 
with the synthetic spectra (red curves) calculated in the si- 
multaneous fit of all hf components to the observed profiles 
(blue histograms). The residuals "observed data - model" 
are depicted beneath each NH 3 spectrum in Fig. B.l. 

The results of the fits to the low-resolution NH3 spectra 
obtained at different positions and the corresponding phys- 
ical parameters are given in Table B.l. The estimate of the 
excitation temperature is made from the (1,1) transition 



assuming a beam filling factor rj = 1 since both the major 
and minor axes of the core exceed the angular resolution. 
At two offsets (0", 0") and (0", -40"), we calculated the ro- 
tational and kinetic temperatures from the two transitions 
(1,1) and (2,2), and from Eq.(A.21) we determine the gas 
densities nn 2 ~ 1-2 x 10 4 and 2.5 x 10 4 cm~ 3 , respectively, 
for a uniform cloud coverage (Col. 10 in Table B.l). 

Assuming spherical geometry, we find a core mass 
M ~ 0.5M Q (the mean molecular weight is 2.8). On the 
other hand, the virial mass [Eq.(A.26)] calculated from the 
linewidth Aw and the core radius r is M v ; r <~ 1.0M Q . The 
observed difference in the masses could be due to devia- 
tions from a uniform gas density distribution and a core 
ellipticity. 

The total ammonia column density of AT tot <~ 1.3 x 10 15 
cm -2 gives the abundance ratio [NH 3 ]/[H 2 ] ~2x 10~ 7 at 
the core center (d <~ 0.08 pc) which is consistent with other 
sources (e.g., Dunham et al. 2011). 

The high resolution NH 3 spectrum (Fig. B.l) can be 
used to estimate the input of the turbulent motions to the 
line broadening. The measured linewidth Av = 0.27 km s _1 
and the thermal width Aw t h ~ 0.15 km s" 1 at Tki n ~ 9 K 
(Table B.l) give er turb ~ 0.09 km s _1 . 

Figure lb shows the NH3(1,1) intensity (grey contours) 
and the radial velocity V LSR (color map) structure in L694- 
2. The radial velocities within the central 115" x 65" re- 
gion do not show significant variation in both magnitude 
and direction, suggesting a simple solid-body rotation. The 
velocity gradient along the major axis (P.A. = 127°) is 
grad a U = (V — (V^ SR ))/0.5a <~ ±1 km s _1 pc -1 , and along 
the minor axis grad b U = (V — (V LSR ))/0.5b ~ and +1 
km s _1 pc -1 . Towards the edges N-W (redshift) and S-E 
(blueshift), ~ 50" off the center, the radial velocity mono- 
tonically changes by ±0.06 km s _1 . Being attributed to 
cloud rotation, such systematic trends in Vl SR correspond 
to an angular velocity cf> w 4 x 10~ 14 s _1 . 

For a self-gravitating rigidly rotating sphere of constant 
density p, mass M, and radius r, the ratio of rotational to 
gravitational potential energy is 



E, 



rot 



u 



3M(j> 2 



l6ir 2 Gp 2 r 3 AirGp 



or 



' = 2.55 x 10~ 3 <£L 14 /n 4 , 



(2) 



(3) 



where 0_i4 is the angular velocity in units of 10 14 s 1 and 
11,4 is the gas density in units of 10 4 cm~ 3 (Menten et al. 
1984). 

Equation (3) and the estimates of the angular velocity 
and the gas density obtained above show that j3 < 0.016 
and, hence, the rotation energy is a negligible fraction of 
the gravitational energy at this stage of evolution of the 
low mass starless molecular core KaOl. 

There is a different measure of the influence of cloud ro- 
tation upon cloud stability through the ratio between the 
rotational and combined thermal and non-thermal (turbu- 
lent) virial terms (e.g., Phillips 1999): 



/3' = 7.04 x 10~ 2 r 2 <j> 2 _ 14 Av~ 



(4) 



where r is the cloud radius in pc, and Av the linewidth 
(FWHM) in km s _1 . The influence of rotation and turbu- 
lence are comparable when the stability parameter (3'(r) w 
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1. For the core KaOl, (3' ss 0.02 which means that turbu- 
lence greatly exceeds the contribution due to rotation in 
determining cloud stability. 

3.2. Dobashi 279 P6 (Do279 P6) 

The source Dobashi 279 P6 (Do279P6, herein) shows a very 
complex structure in the NH3 emission (Fig. 2a) with many 
small clumps located along the main axis (P.A.= 56°, a « 
9', or ~ 0.5 pc at D ~ 203 pc). We list the coordinates 
of five bright NH 3 peaks, the peak temperatures T MB , their 
radial velocities V LSR , and the linewidths Av in Table 1. The 
peaks are labeled from a to e in Fig. 2b and in Table 1 in 
accord with decreasing main beam brightness temperature, 

The location of Do279P6 coincides with a star formation 
region in the Serpens cloud, which is known under the name 
Cluster B (Harvey et al. 2006), or Scr G3-G6 (Cohen & 
Kuhi 1979), or Serpens NH 3 (Djupvik et al. 2006). The area 
around Ser G3-G6 was mapped in the NH 3 (1,1) emission 
line by Clark (1991) who found two ammonia condensations 
on each side of the complex, Ser G3-G6NE and Ser G3- 
G6SW (the map is not published; observations obtained 
with the Effelsberg 100-m telescope). The coordinates of 
these condensations coincide, respectively, with the a and 
7 ammonia peaks, whereas the close group of G3-G6 stars 
lies in the field of the S peak. 

The Do279P6 cloud with identified YSO candidates is 
shown in Fig. 2c in which open circles of different colors 
mark Class I (red), "flat" (yellow) 2 , and Class II (blue) 
candidates. Class II sources are more evolved and older 
than Class I and, thus, the NH3 core might be already 
more dispersed. So, what we observe at the S — W off- 
set 7 = (—160", —80") is consistent with this expectation - 
only "flat" and Class II sources are located therein. On 
the other hand, all Class I sources (ID numbers B13 and 
B15-B18 in accord with Table 4 in Harvey et al. 2006) are 
detected in the vicinity of the core center a = (0", 0") which 
has the highest gas density, nn 2 ~ 2.3 x 10 4 cm~ 3 , and the 
maximum intensity of the ammonia emission T MB = 2.9 K 
(Table B.2). We note that the source B18 lies within 8" off 
the NH3 peak, and B13 - the most distant - at <~ 80" along 
the N-W direction. It is also noticeable that the angular size 
of the central condensation is two times smaller then that of 
the S-W clump, i.e., the region is probably less dispersed. 
The two other condensations with "flat" and Class II type 
sources are located at the offsets 5 = (—120",— 40") and 
(3 = (40", 40") with, respectively, five and one embedded 
YSO(s). The only starless NH 3 region detected at the e 
condensation e — (80", 80") has the narrowest lincwidth 
AV ~ 0.4 km s _1 (Table 1) and the lowest kinetic tem- 
perature Tkin — 10 K (Table B.2) which are comparable to 
KaOl — the isolated dense core (Sect. 3.1). 

Figure B.2 shows the high-resolution ammonia spec- 
tra (channel spacing 0.015 km s" 1 ) obtained at three NH 3 
peaks: a, (3, and 7. The corresponding linewidths are the 
following: Av ss 0.75, 0.63, and 0.96 km s _1 . Given the de- 
rived kinetic temperatures (Table B.2), we can calculate the 



2 A "flat" category is between Classes I and II (Andre & 
Montmerle 1994; Greene et al. 1994; Deharveng et al. 2012). 



thermal contribution to the linewidth and the non-thermal 
velocity component: 

Ai = A„*- 8 l»2^, (5) 

where (fcTki n /17TOH~) is the thermal broadening due to NH 3 , 
and mH is the mass of the hydrogen atom. 

The calculated non-thermal velocity dispersion, cr NT = 
Aw NT /(2v / 2 In 2), of about 0.3 km s _1 is comparable to the 
large scale motions observed in the cloud. Namely, the N — 
E part is moving towards us (with respect to the systemic 
velocity of the cloud, (V^ SR ) = 7.986 km s _1 ) with a radial 
velocity V^ SR w —0.4 km s^ 1 and the S — W part in the 
opposite direction with similar radial velocity (Fig. 2b). In 
the middle of the cloud there is a narrow zone along the 
cut Aa = —50" with positive Vl SR rs 0.3 km s _1 , which is 
bracketed by two zones with negative V^ SR w —0.3 km s _1 . 
The whole structure may be due to the differential rotation 
of different clumps distributed along the main axis of the 
cloud. 

The central ammonia condensation consists of two 
clumps a and (3 separated by ~ 60". The position angle 
of the major axis of the brighter a clump is P. A. ~ 50°. 
The deconvolved diameters of the major and minor axes 
are a ~ 130" and b ~ 60" (FWHP). The angular sizes of 
the j3 clump are less certain since it is not completely re- 
solved. The physical parameters were estimated assuming 
a filling factor 77 = 1 for both the a and j3 clumps and thus 
T cx , A^(NH 3 ), and n(H2), listed in Table B.2, correspond to 
their minimum values (for details, see Appendix A). With 
the mean diameter of ~ 90" (or 0.09 pc at D = 203 pc) the 
mass of this double clump is M a ^ > 0.5M©, and the am- 
monia abundance is [NH 3 ]/[H2] ~2x 10~ 7 . The estimate 
through virial masses is less certain in this case since YSOs 
are embedded within the clumps. 

Another double peak is formed by the 7 and (5 cores. 
The major axis of the brighter 7 condensation (P.A. = 0°) is 
113" and the minor axis is 69" (deconvolved, FWHP). The 

5 clump is not completely resolved which gives, as for the 
previous pair, the lower bound on the mass M 7; ^ > 0.4M Q 
and the abundance ratio [NH 3 ]/[H 2 ] — 1 x 10~ 7 . 

The weakest peak, e, separated by <~ 60" from f3, is a 
completely unresolved ammonia clump with an angular size 

6 < 40". It shows the largest optical depth, m ~ 11, and 
the lowest kinetic temperature, Tk; n ~ 10 K (Tabic B.2). 
Other physical parameters listed in Table B.2 were calcu- 
lated at rj = 1, i.e., they are the minimum values. The 
filling factor 77 is restricted in this case between 0.2 and 1 
(see Appendix A), which provides the following boundaries: 
4 < T cx < 10 K, 1.3 x 10 15 < iV(NH 3 ) < 3.2 x 10 15 cm" 2 , 
0.8 x 10 4 < n(H 2 ) < 4.9 x 10 4 cm" 3 . The mass of the e 
clump is less than O.O6M if n(H 2 ) — 3 x 10 4 cm" 3 and its 
diameter d < 0.04 pc. 

Of course, the derivation of relative ammonia abun- 
dances and masses is not very accurate and based on a 
number of model assumptions like spherical geometry, ho- 
mogeneous gas density, the distance to the target of 203 pc, 
etc. Nevertheless, the substellar mass of the e clump may 
lie in the range of the lowest possible masses with which a 
brown dwarf can form: 0.03 — 0.08M Q (e.g., Oliveira et al. 
2009). 
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3.3. Dobashi 279 P8 (Do279 P8) 

The molecular cloud Dobashi 279 P8 (Do279P8, herein) 
is an example of a starless core with no embedded IRAS 
point sources or a pre-main-sequence star. Our mapping 
reveals a slightly elongated structure of this core in the 
N — S direction (Fig. 3) with an FWHP angular size of 
<~ 45" x 57" (deconvolved), or ~ 0.04 x 0.06 pc 2 linear size 
at D = 203 pc. The coordinates of the NH 3 peak, the peak 
temperature T MB , its radial velocity V LSR , and the linewidth 
Av (FWHM) are listed in Table 1. The NH 3 spectrum is 
shown in Fig. B.3. We did not detect any emission at the 
expected position in the (2,2) transition and, therefore, only 
an upper limit on the rotational temperature is given in 
Table B.l. The la upper limit on T rot is very low, < 10 K. 
A noticeable feature of the Do279P8 spectrum is a very 
narrow linewidth of the hfs transitions, Av ~ 0.26 km s _1 . 
Another remarkable characteristic is slow rotation of the 
cloud illustrated in Fig. 3b: at a radius of <~ 20" off the 
center (« 0.02 pc), the tangential velocity does not exceed 
±0.04 km s _1 , which leads to an angular velocity cb <~ 6 x 
10~ 14 s _1 . The axis of the velocity centroid is slightly tilted 
with respect to the main N — S axis of the NH 3 (1,1) map. 

The measured excitation temperatures at the offsets 
(0", 0") and (0", -40") are close to the upper limits on the 
rotational temperatures (see Table B.l). This may imply 
that the (1,1) transition is almost thermalized and ammo- 
nia emission arises from an extended structure which fills 
the telescope beam homogeneously, r\ ss 1. At our angular 
resolution of 40", the NH 3 (1,1) hfs components do not show 
any blueward or redward asymmetry indicating contract- 
ing or expanding motions, respectively. Most probably, the 
core Do279P8 is in an early stage of evolution. According 
to Lee & Meyers (2011), the starless cores evolve from the 
static to the expanding and/or oscillating stages, and fi- 
nally to the contracting cores. The suggestion that the core 
is in an initial stage of evolution is also supported by low 
values of (3 < 0.1 and (3' s=s 0.01. In this estimation, we 
use n gas > 1 x 10 4 cm" 3 ([NH 3 ]/[H 2 ] - 5 x 10~ 7 ) as a 
lower bound on the gas density. With the observed diameter 
d <~ 50" (deconvolved) and assuming spherical geometry, 
the lower limit on the mass of Do279P8 is M ~ O.O6M . 

3.4. Dobashi 279 P12 (Do279 P12) 

The source Dobashi 279 P12 (Do279P12, herein) is a dark 
cloud with the opacity class 4-5 in accord with Lynds 
(1962). The cloud is known under the names L572 (Lynds 
1962), BDN 31.57+5.37 (Bernes 1977), or the Serpens dark 
cloud which is a site of active star formation (Strom et al. 
1974). A large number of articles has been devoted to study- 
ing this cloud in all spectral ranges from radio to X-ray (for 
a review, see, e.g., Eiroa et al. 2008). 

Single-dish ammonia observations of the Serpens dark 
cloud were started by Little et al. (1980) and by Ho & 
Barrett (1980). Later on, Ungerechts & Glisten (1984) 
mapped an extended ammonia structure of this cloud with 
the Effelsberg 100-m telescope. Ammonia interferometric 
observations with the VLA by Torrefies et al. (1992) were 
directed to the study of the embedded high-density molec- 
ular gas around the triple radio source Serpens FIRS 1 - 
a typical Class protostar (Snell & Bally 1986; Rodrigues 
et al. 1980). Further interferometric (VLA) and additional 
single dish (Haystack 37-m, Effelsberg 100-m) observations 



of FIRS 1 by Curiel et al. (1996) revealed high-velocity 
ammonia emission (up to <~ 30 — 40 km s _1 from the line 
center) aligned with the radio continuum jet. 

The results of our Effelsberg observations are shown in 
Figs. 4-6. The NH 3 (1,1) map (Fig. 4) exhibits a complex 
structure with three bright peaks, T MB ~3K, located NW- 
SE at the offsets a = (160", -40"), fi = (200", -120"), and 
7 = (40", 40") (see Table 1 and Fig. 5, panel B). The a 
peak lies very close to the maximum of the NH 3 emission 
detected by Little et al. (1980), whereas the (3 and 7 peaks 
are spatially coincident with positions of the two peaks lo- 
calized by Ho & Barrett (1980). The measured NH 3 profiles 
show in this case a double structure in both the (1,1) and 
(2,2) transitions at four positions (Fig. B.4, Table B.3) and 
a split (1,1) line at another three positions where the (2,2) 
line was not detected. This multicomponent structure has 
not been resolved previously in the ammonia observations 
cited above. Since the splitting is also seen in the optically 
thin ammonia lines, the observed asymmetry of the NH 3 
profiles has clearly a kinematic origin. 

We used a two-component model, Eq. (A. 7), to fit the 
data. The result is illustrated in Fig. B.4 by the red lines, 
the detected velocity subcomponents are labeled by letters 
A and B. The seven offsets with the split NH 3 lines are 
marked by white crosses on the velocity maps shown in 
Fig. 5. The cluster of the white crosses at the SE part 
of Do279P12 outlines a small area of the bipolar outflow 
with an angular extent of ~ 40" x 80". Table B.3 displays 
the physical parameters estimated at different offsets. The 
mean radial velocities of the components A and B averaged 
over the seven offsets are (Va) = 7.27 ±0.11 km s _1 and 
(Vb) = 8.54 ± 0.08 km s , and the velocity difference is 
AVb_ a = 1.3 ± 0.1 km s . Panel A in Fig. 5 shows that 
the gas surrounding the SE double component region has 
predominantly a negative radial velocity with respect to 
(Ksr), whereas the gas shown in panel B has a positive 
Ksr - The observed velocities converge to the mean (Ksr) 
at the cloud's outskirts. 

Figure 5 shows also a velocity gradient along the E-W 
direction around the isolated double component region at 
(80", 0") with V A = 7.89 km s" 1 and Vb = 8.92 km s" 1 . 
The adjoint regions at (120", 0") and (40", 0") have inter- 
mediate velocities of Vl SR = 8.14 and 8.42 km s _1 , re- 
spectively, which become closer to the mean radial veloc- 
ity of the cloud (K SR ) = 8.099 km s" 1 at larger offsets: 
F LSR = 7.90 and 8.04 km s' 1 at (280", 0") and (-40", 0"), 
respectively. 

The measurement of the excitation temperatures reveals 
T ex essentially lower than T rot at all 23 offsets listed in 
Table B.3. As discussed above, this may indicate dumpiness 
of the ammonia distribution and 77 < 1. In Fig. 4b, we also 
show the Herbig-Haro (HH) objects and far-infrared con- 
tinuum sources taken from Davis et al. (1999), the submm- 
continuum sources detected by SCUBA (Di Francesco et al. 
2008), and the position of the triple radio source FIRS 1 
(Snell & Bally 1986) which coincides with the brightest in- 
frared object SMM1. The closest position of the split NH 3 
lines to the bipolar outflow in FIRS 1 is (80", 0") which cor- 
responds to a projected angular distance of 6 = 45" from 
the SE knot of the radio source FIRS 1 (see also Fig. 13 
in Snell & Bally 1986). All other splittings are located in 
the vicinity of the /3 peak, around the SMM2 and SMM11 
sources. 
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The observed linewidths Av of the split NH3 compo- 
nents are in the range between 0.62 km s _1 and 0.90 km s _1 . 
However, larger linewidths (Av > 1.0 km s _1 ) measured at 
different offsets (Table B.3) may imply unresolved subcom- 
ponents of the ammonia lines due to the moderate spectral 
resolution of our database. Gas flows and jets usually ac- 
company the mass accretion processes in YSOs and proto- 
stars. From this point of view the split and/or widened line 
profiles of NH3, — a tracer of the dense cores, — may indi- 
cate the first stages of the accretion processes in protostars, 
where jets of HH objects have not yet been formed. 

It may be illustrative to compare our Vlsr values with 
those measured with the same telescope and beam size 
in 1977-1981 (Ungcrechts & Giistcn 1984). In total, we 
have 29 overlapping positions shown in Fig. 6. The cal- 
culated differences AV = Vigsi — V2011 arc marked by 
dots with 3(J error bars. The red points indicate devia- 
tions from zero which are larger than 3ct. At two positions, 
(Aa, A<5) = (80", 0") and (160", -120"), we detected dou- 
ble NH 3 profiles. Among these 29 positions, significant de- 
viations from zero are revealed at 6 offsets, but other points 
show no variations over this time scale. The ammonia line 
splitting was not seen at (80", 0") in 1981, although the 
linewidth was Auiggi = 0.85±0.20 km s _1 , which is compa- 
rable to that given in Table B.3. However, at (160", -120") 
the linewidth was Awiggi = 1.68 ± 0.17 km s _1 , i.e., two 
times larger than the value from Table B.3. We note that 
the weighted mean of the linewidth at the overlapping po- 
sitions is (Aw) i98i = 0.93 km s _1 , and the dispersion of the 
sample is cr& v — 0.18 km s _1 . The enhanced Av(NH3) at 
(160", —120") observed 30 years ago was probably caused 
by the unresolved structure of the ammonia lines. The max- 
imum deviation AF LSR <~ 1 km s _1 (Fig. 6) constrains the 
gas acceleration in an outflow at the level of V < 0.03 
km s _1 yr -1 . This means that the high-velocity ammonia 
emission up to ~ 30 — 40 km s _1 from the line center, ob- 
served, for example, in FIRS 1 (Curiel et al. 1996), could be 
reached in a period of <~ 1000 yr. This is comparable to the 
typical time scale of 10 3 yr — a scale for the episodic ejec- 
tion of material from protostars (e.g., Ioannidis & Froebrich 
2012). The mass of the protostar should be < O.3M in 
order to eject the gas clump at the 40 AU radius with the 
acceleration V ~ 0.03 km s _1 yr -1 . 

3.5. Dobashi 321 P2 (Do321 P2) 

The source Dobashi 321 P2 (Do321P2, herein) is a fil- 
amentary dark cloud with the opacity class 6 in accord 
with Lynds (1962). The cloud is known under the names 
L673 (Lynds 1962), BDN 46.22-1.34 (Bernes 1977), or P77 
(Parker 1988). Do321P2 was observed in ammonia lines by 
Anglada et al. (1997) and Tobin et al. (2011). 

Our NH3 observations are shown in Fig. 7. The ammo- 
nia structure consists of five main subcondensations. Their 
reference positions are listed in Table 1. The brightest a and 
(3 peaks with T MB ~3K are located near the sources IRAS 
19180+1114 (SMM2) and IRAS 19180+1116 (SMM1), re- 
spectively. The ammonia peak detected by Anglada et al. 
has an offset (—123", 83") and is resolved into two clumps 
a and (3 in our observations (see Fig. 7). We also detected 
three weaker ammonia clumps 7, 6, and s with T MB <~ 
1.7, 1.6, and 1.3 K, respectively. The 7 and e condensations, 
localized within the confines of the older map of Anglada et 



al. (1997), peak near the source IRAS 19181+1112, whereas 
the S clump was not resolved by Anglada et al. (1977) be- 
cause of a too coarse angular resolution. 

Interferometric observations of a small region around 
SMM2 (~ 100" x 100") in lines of N 2 H+ and NH 3 (Tobin 
et al. 2011; see Fig. 14) show a substantial amount of ex- 
tended ammonia emission which closely traces the spacial 
distribution of N2H" 1 ". The measured NH3 radial velocity 
V LSR = 6.93 km s _1 , the linewidth Av — 0.50 km s _1 , and 
the column density iV(NH 3 ) = 7.2 x 10 14 cm" 2 (Table 8 
in Tobin et al.) are in a good agreement with our values 
for the a peak (Table 1): V LSR = 6.88 km s" 1 , Ai; = 0.49 
km s _1 , and NQSIHs) = 6x 10 14 cm~ 2 . However, the es- 
timates of the excitation temperature at this position con- 
tradict sharply: T ex = 33.3 K (Tobin et al.) and 6.8 K (our 
value from Table B.l). We detected NH 3 (1,1) and (2,2) at 
the a and j3 peaks (Fig. B.5). Their profiles are well de- 
scribed by a single-component model (red lines in Fig. B.5, 
model parameters are given in Table B.l). Both spectra 
show the same T ex , which is also consistent with the values 
measured in other dense molecular cores (see Tables B.2- 
B.4). 

The NH 3 velocity map (Fig. 7b) shows irregular struc- 
ture with distinct pockets of redshifted and blucshiftcd am- 
monia emission around the clusters of YSOs at the positions 
of SMM1 and SMM2 (Fig. 7c). The observed kinematics are 
caused by large-scale motions in this star-forming region 
which cannot be interpreted unambiguously. 

The apparent diameters (beam deconvolved) of the two 
NH3 condensations d a — 60" and dp — 40" and the gas 
densities listed in Table B.l allow us to estimate masses of 
these clumps assuming spherical geometry: M a ~ O.3M0, 
and Mp ~ O.IMq. The relative ammonia abundance in 
both clumps is [NH 3 ]/[H 2 ] - 2 x 10~ 7 . We note that the 
corresponding virial masses of these clumps are consider- 
ably larger (M^ ir ~ 2M , and M? r ~ 1M ) and probably 
affected by the masses of YSOs embedded in the cloud. 

3.6. Kawamura 05 (Ka05) 

The source Kawamura 05 (Ka05, herein) is a Bok globule of 
the opacity class 6 (Lynds 1962). The cloud is known under 
the names B335 (Barnard 1927), L663 (Lynds 1962), or 
CB199 (Clemens & Barvainis 1988). Being the prototype of 
an isolated, star-forming dark cloud, B335 has been studied 
at many wavelengths. 

Single-pointing ammonia observations towards B335 
were carried out by Ho et al. (1977, 1978), Myers & Benson 
(1983), Benson & Myers (1989), and Levshakov et al. 
(2010). The first mapping of this core in the NHs(l,l) and 
(2,2) lines with the Effelsberg telescope revealed an elon- 
gated N-S structure (Menten et al. 1984, M84 herein). The 
kinetic temperature measured at four positions, where the 
NH 3 (2,2) line was observed, is uniformly distributed over 
the core region, Tkin = 10-12 K. 

We mapped the core Ka05 at nine offsets around 
the reference point a(J2000) = 19:37:01.3, <S(J2000) = 
+07:34:29.6 which is shifted off the center position in M84 
at (Aa, A6) = (-24.0", -42.8"). The results of our ob- 
servations are shown in Fig. 8. Two positions (0",0") and 
(0",-40") with T MB = 2.5 K and 2.1 K (Table 1) coin- 
cide, respectively, with the offsets (20", 40") with T MB = 
3.4 K and (20", 0") with T MB = 2.8 K in M84. In between 
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these positions, at (20", 20"), the brightness temperature 
was even higher, T MB — 4.1 K (M84), but we did not ob- 
serve this point with our larger step size of 40", as well 
as the point (20", 60") where T MB = 3.1 K. However, we 
measured NH 3 (1,1) at (0",40") with T MB = 0.6 K which 
corresponds to the point (20", 80") not listed in Table 1 in 
M84. The opposite marginal point (20", -20") along the 
cut Aa = 20" in M84 has T MB = 0.8 K. This indicates 
sharp drops of the brightness temperature at the edges of 
the disk-like envelope surrounding the protostar. 

The spectra with the measured (1,1) and (2,2) lines 
are shown in Fig. B.6. The radial velocities at two off- 
sets (0",0") and (0",— 40") are the same as they were in 
1982/83 (cf. Fig. 6, Sect. 3.4, where velocity shifts at some 
positions in the map of Do279P12 were detected for the 
same three decades). The linewidths (Table B.l) are also 
consistent with M84. At these offsets, the measured ki- 
netic temperatures, excitation temperatures, total optical 
depths, and column densities are in line with M84. 

The NH3(1,1) velocity map is shown in Fig. 8b. The 
velocity distribution resembles a rigid body rotation with 
the tangential velocity increasing from w —0.040 km s _1 
(NE) to w +0.040 km s" 1 (SW) at cloud flanks. The radius 
of the NH 3 disk-like envelope is about 40" in accord with 
our and the M84 maps. The distance to B335 was recently 
redetermined to be only D = 90 — 120 pc (Olofsson & 
Olofsson 2009). This gives us a radius r ~ 0.02 pc, and 
an angular velocity <p ~ (6 - 8) x 10~ 15 s" 1 ((3 < 0.001), 

which is close to the value of M84, <j) ~ 1.5 x 10~ 14 s _1 . 
The corresponding stability parameter j3' « 2 x 10~ 4 . With 
the apparent diameter d <~ 0.04 pc, the relative abundance 
ratio is [NH 3 ]/[H 2 ] - 4 x 10~ 7 . 

To compare with previous estimates of the mass of the 
gas giving rise to ammonia emission, we assume spheri- 
cal geometry and for the gas density riH 2 ~ 2 x 10 4 cm~ 3 
(Table B.l) to obtain M ~ O.O5M — 40 times lower than 
the value reported by Benson & Myers (1989) who used a 
core radius R = 0.11 pc (D = 250 pc) and a gas density 
riH 2 ~ 6 x 10 3 cm~ 3 . 

If the centrifugal forces prevent the disk-like structure 
from collapse, then the mass of the center source is given, 
within a factor of 2, by (e.g., Frerking & Langer 1982): 

M c = V t 2 R/G, (6) 

where R is the distance of the material from the source and 
G is the gravitational constant. For V t = 0.040 km s _1 , it 
gives M c <~ 0.03M Q . This estimate indicates that Ka05 is 
a low-mass protostellar source which agrees with the mea- 
surement of the central source mass M c ~ 0.04M Q from 
SMA observations of the C 18 0(2-l) emission (Yen et al. 
2010). 

3.7. Dobashi 279 P7 (Do279 P7) 

The dense molecular cloud Dobashi 279 P7 (Do279P7, 
herein) is the next example of a starless core from our 
list. This is a newly discovered dense clump which is so far 
poorly studied. No embedded IRAS sources or pre-main- 
sequence stars are known within 80" off the core center. 
The NH 3 (1,1) intensity distribution shown in Fig. 9 has 
approximately circular symmetry. The half-power dimen- 
sion of the core is 80" x 80", or, correcting for our beam, 



70" x 70", i.e., the source is resolved in both right ascen- 
sion and declination. Assuming a distance to the source of 
D w 203 pc (see Sect. 1), its radius is about 0.03 pc. The 
parameters for the density peak at (—80", 40") are given 
in Table 1. The NH 3 spectra are shown in Fig. B.7. We 
detected the (2.2) transition at four positions which allow 
us to estimate T rot = 11.5 - 13.5 K and T kin = 12.4 - 14.9 
(Table B.l). These values exceed slightly the temperatures 
measured in the starless core KaOl, but the linewidths in 
Do279P7, are significantly larger than the linewidths in 
KaOl and Do279P8 - another starless core but without de- 
tected (2,2) transition (see Table B.l). The wider lines in 
Do279P7 indicate a higher level of turbulence in this core 
with the ratio of turbulent to thermal velocity dispersion 
of A^turb/Awth ~ 5. The profiles of the NH 3 (1,1) hfs com- 
ponents do not show any blueward or redward asymmetry 
indicating contracting or expanding motions. However, the 
NH 3 intensity map shows a filament-like structure in the S- 
E direction (P.A. w 135°) with V LSR ~ 0.8 km s" 1 (Fig. 9). 

The core exhibits a linear velocity gradient from the 
southern (V LSR ~ 0.0 km s _1 ) to the northern edge (Ksr 
~ -0.8 km s" 1 ) that can be interpreted as rigid body rota- 
tion around the axis Ad w 20", P.A. w 90° (Fig. 9b). The 
measured tangential velocity of ±0.4 km s _1 atrw 0.03 pc 
leads to an angular velocity of w 5 x 10~ 13 s _1 . With 
rtH 2 ~ 1 x 10 4 cm~ 3 (Table B.l) and assuming spherical 
geometry, one obtains a core mass of M <~ O.2M . The 
estimates of the angular velocity and the gas density give 
j3 <~ 6 which shows that the rotation energy is larger or 
comparable to the gravitational energy of this starless core. 
The parameter /?' w 0.3 is also rather high as compared 
with other isolated cores KaOl, Do279P8, and Ka05 where 
/?' <C 1. The measured column density N <~ 6 x 10 14 cm~ 2 
and the linear size d <~ 0.06 pc indicate a relative ammonia 
abundance of [NH 3 ]/[H 2 ] - 3 x 10~ 7 . 

3.8. Serpens South 3 (SS3) 

The source Serpens South 3 (SS3, herein) is a dense core 
located in an active star-forming region discovered by 
Gutermuth et al. (2008) in Spitzer IRAC mid-infrared imag- 
ing of the Serpens- Aquila rift. The central part of this re- 
gion, near a = 18 h 30 m 03 s , S = -02° 01' 58.2" (J2000), 
which is called Serpens South, is composed of 37 Class I 
and 11 Class II sources at high mean surface density (> 430 
pc~ 2 ) and short median nearest neighbor spacing (0.02 pc). 
The overall structure of the star-forming region exhibits 
a typical filamentary dense gas distribution elongated for 
about 14' in the NW-SE direction. 

For the first time, we mapped the core SS3 in the 
NH 3 (1,1) and (2,2) lines around the reference point a = 
18 h 29 m 57.1 s , S = -02° 00' 10.0" (J2000) in the area 
Aa x AS = 160" x 360" (Fig. 10). The nominal center 
of the NH 3 map is shifted by 2'29" (P.A. w -40°) with 
respect to the center of Box 1 in the CO map of Nakamura 
et al. (2011, herein Nil) and coincides with the intensity 
peak of the CO outflow lobe marked by R2 in Nil. The 
angular size of the CO lobe within the 15 K km s _1 con- 
tour is about 70" x 60" (Figs. 9 and 15 in Nil), and the 
whole region lies in the lowest NH 3 intensity 'valley' run- 
ning in E-W direction along the cut AS = 0" in Fig. 10. 
To the south of this valley, at AS w —140", an intensity 
peak of NH 3 (a peak in Table 1 and Fig. 10b) occurs at 
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the northern part of another CO outflow lobe marked R3 
in Nil. These two redshifted lobes R2 and R3 appear to 
have the blucshifted counterparts Bl and B4 originating 
from a compact 3 mm source — a deeply embedded, ex- 
tremely young Class protostar (Nil). It is interesting to 
note that the configuration of the NH3 valley just follows 
the lowest dust emission distribution on the 1.1 mm con- 
tinuum image of this area shown in Fig. 15(b) in Nil. To 
the north and to the south of the valley there are regions 
of relatively strong 1.1 mm emission which are correlated 
with strong NH 3 emission. The three NH 3 peaks a, (3, and 
7 listed in Table 1 and labeled in Fig. 10b have very high 
brightness temperatures T MB = 5.0, 3.1, and 3.0 K. 

Examples of the observed NH3(1,1) and (2,2) spectra 
are shown in Fig. B.8. Here again, as was the case of the core 
Do279P12, we find split NH 3 profiles. The apparent asym- 
metry of the ammonia lines is observed in optically thin hf 
components and, thus, cannot be due to radiative trans- 
fer effects (self-absorption). The components A and B, de- 
tected at two offsets (see Table B.4), have the radial veloc- 
ities (V A ) = 7.865 ±0.005 km s" 1 and (V B ) = 6.856 ±0.010 
km s _1 . The map of the velocity field is shown in Fig. 10b 
where the positions of the split NH3 profiles are labeled by 
the white crosses. The closest Class protostar to these 
positions is the MAMBO mm source MM10 (Maury et al. 
2011) marked by the yellow circle in Fig. 10c at the offset 
(38", —53"). However, since the whole region is very com- 
plex, it is not clear which source drives this presumably 
bipolar molecular outflow causing the observed splitting of 
the NH 3 profiles. 

The measured linewidths listed in Table B.4 show that 
the values At; > 0.8 km s _1 are symmetrically distributed 
around the positions of the split NH 3 profiles in the E-W 
strip of about 120" width, cut by the vertical edges of our 
map. The strip of enhanced linewidths (see Fig. B.12) may 
be associated with the E-W bipolar outflow, thus giving 
rise to the local redshifted and blueshifted patches of the 
velocity map shown in Fig. 10b. 

It is seen from Table B.4 that the values of the mea- 
sured parameters fluctuate considerably from one position 
to another. Namely, the excitation temperature ranges from 
3.3 K to 8.4 K, T rot — from 8.6 K to 14.8 K, T kin — from 
8.9 K to 16.8 K, the gas density nn 2 — from 0.3 x 10 4 cm~ 3 
to 4.2 x 10 4 cm~ 3 , the total optical depth m — from 1.5 
to 20.0, and the column density JV(NH 3 ) — from 0.2 x 10 15 
cm~ 2 to 2.8 x 10 15 cm~ 2 . We note that the highest gas den- 
sity of «h 2 ~ 4.2 x 10 4 cm -3 was measured at two positions 
where we observed the split NH 3 lines profiles. The largest 
optical depths with 13.2 < m < 20.0 are localized along 
the N-S cut Aa = 0" between AS = 80" and 200" which 
passes through the 7 and f3 peaks in the region where three 
starless core candidates (red circles in Fig. 10c) were de- 
tected by Maury et al. (2011). The maximum (16.6 K and 
16.8 K) and minimum (8.9 K) kinetic temperatures are ob- 
served along the N-S cut at Aa = 80". The minimum of 
the kinetic temperature correlates with the minimum of the 
brightness temperature, T MB = 0.5 K, whereas T MB = 1.2 K 
at Tkm= 16.8 K (the foot of the a peak) and T MB = 5.0 K 
at Tki n = 16.6 K (the top of the a peak). This gradient of 
?kin (8.9 — 16.8 K) is probably due to the presence of two 
YSOs localized at the top and foot of the a peak shown in 
Fig. 10b. 



3.9. Dobashi 243 P2, 279 P18, 279 P13, and 321 PI 

The sources Do243P2, Do279P18, Do243P13, and Do321Pl 
(Table 1) were observed at their central positions and at 40" 
offsets toward the four cardinal directions. Only a weak 
NH 3 (1,1) emission (T MB < 1 K) is seen toward each of 
them. The corresponding spectra are depicted in Fig. B.9. 
For Do321Pl, two possible detections of the ammonia emis- 
sion are shown at the bottom panels, but the noise is 
high and the identification of the NH 3 (1,1) line is not cer- 
tain. The measured radial velocity of V LSR = 13.5 km s _1 
for Do321Pl is somewhat larger than that measured for 
the other sources from Table 1 for which Vl SR ranges be- 
tween 6.0 km s _1 and 9.6 km s _1 , however, the 12 CO(T 
0) 115.3 GHz, 13 CO(1-0) 110.2 GHz, and C 18 O(l-0) 109.8 
GHz emission lines were also detected at V LSR ~ 13 km s _1 
with the Dclingha 14-m telescope (Wang 2012). 

No submillimeter-continuum sources have been found 
within 2' off the core centers of Do279P18, Do243P13, and 
Do321Pl, but a SCUBA point-like source J180447.8-043208 
(Di Francesco et al. 2008) lies in the field of Do243P2 at the 
angular distance 9 = 41", P.A. = —30° from the nominal 
center (0",0"). 

3.10. Sources without detected A//-/ 3 emission 

A list of 37 targets from our dataset with undetected NH 3 
emission is presented in Table 2. Columns 3 and 7 give the 
noise level per channel width of 0.077 km s _1 on a T MB 
scale and the mapped area, respectively. In the fields of 
some targets there were found the following objects. 

Do279P21: a point-like infrared source IRAS 18191-0310 
(Cutri et al. 2003) at 6 = 80", P.A. = -174°. 

Do279P17: a point-like infrared source IRAS 
18351+0010 (Cutri et al. 2003) at 6 = 17", P.A. = 
74°. 

Do296P3: a point-like submillimeter-continuum source 
J185121. 1-041714 (Di Francesco et al. 2008) at 6 = 29", 
P.A. = -173°. 

Kal5: a galaxy gJ194852.5-103056 (Skrutskie et al. 
2006) at 6 = 64", P.A. = 99°. 

Kal7: a point-like infrared IRAS 19470-0630 (Porras et 
al. 2003) at 6 = 57", P.A. = 162°. 



4. Summary 

We have used the Effelsberg 100-m telescope to observe the 
NH 3 (1,1) and (2,2) spectral lines in high density molecu- 
lar cores. The targets were preliminarily selected from the 
CO and 13 CO survey of the Aquila rift cloud complex car- 
ried out with the Delingha 14-m telescope. We measured 
a minimum grid of 5 positions with 40" spacing centered 
upon the (0",0") position determined from the maximum 
intensity of the CO maps. We mapped a larger area for 
sources where ammonia emission was detected to delineate 
the spatial distribution of NH 3 . 

In total, we have mapped the first 49 CO sources (from 
<~ 150 targets) and detected 12 NH 3 emitters what gives 
an overall detection rate of ~ 24 per cent which is 3 times 
lower as compared to the Pipe Nebula which was observed 
with the GBT 100-m telescope (Rathborne et al. 2008) 
and which is not too far from the clouds studied here (see 
Sect. 1). The NH 3 sources in our sample represent diverse 
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populations of molecular clouds from isolated and homo- 
geneous starless cores with suppressed turbulence to star 
formation regions with complex intrinsic motion and gas 
density fluctuations. 

The starless cores KaOl, Do279P8, Do279P7, 
Do279P18, Do243P13, and a candidate - - the clump 
e Do279P6 — have radial velocities Vl SR ranging from 
6.0 km s _1 to 9.6 km s _1 with the mean (Vl SR ) = 8.1 
km s _1 , and dispersion a = 1.2 km s _1 . The only starless 
core with beyond the 3a interval around the mean is 
Do321Pl for which K LSR = 13.5 km s _1 . Three cores - 
KaOl, Do279P8, and the clump e Do279P6 — demonstrate 
rather narrow linewidths, Av < 0.4 km s -1 , whereas the 
others have Av <~ 0.5 — 1 km s _1 . The masses of the 
starless cores do not exceed a solar mass and are typically 
(0.1 - O.5)M , except for the cores e Do279P6, Do279P8, 
and Ka05 for which M < 0.1M Q . The measured ammonia 
abundance, X = [NH 3 ]/[H 2 ], ranges between 1 x 10~ 7 and 
5 x 10~ 7 with the mean (X) = (2.7 ± 0.6) x 10~ 7 . 

For three cores — KaOl, Do279P8, and Ka05 — we 
measured low angular velocities (parameters /3 and /?' <C 1) 
and their rotation energies are less than the gravitational 
energy, but Do279P7 — a fast rotator — yields a rotational 
energy comparable to the gravitational energy (/3 « 6 and 
/?' s=s 0.3). It also has a high ratio of the turbulent to thermal 
velocity. 

Two isolated star-forming cores, Ka05 and Do243P2, 
were observed. The cores' radial velocities are Vlsr — 8.4 
km s _1 and 7.0 km s _1 , respectively, and, thus, both of 
them belong to the same group of cores with V LSR <~ 8 
km s . For the former core, we estimate a sub-solar mass 
and thus confirm that Ka05 is a low-mass protostcllar 
source. Both cores exhibit narrow linewidths, Av ~ 0.3 
km s _1 . 

Most intense ammonia emission was observed towards 
four clouds with complex gas density and velocity struc- 
tures which harbor numerous YSOs: Do279P6, Do279P12, 
Do321P2, and SS3. Each of them consists of a number of 
clumps with different angular sizes. The (Vlsr) values show 
that the complex clouds have the same radial velocities 
as the starless cores and isolated star-forming clouds from 
our dataset. It is only one cloud — Do279P12 — where the 
dispersion of bulk motions is comparable to the ammonia 
linewidths. The internal bulk motions in the other three 
clouds are less pronounced and show lower dispersions as 
compared to the NH 3 linewidths. 

Kinematically split NH3 profiles were detected in 
Do279P12 and SS3. The velocity splitting is ~ 1 km s _1 . 
In the former object, the splitting is most probably caused 
by bipolar molecular outflows observed in other molecular 
lines. This object has been observed in ammonia lines with 
the Effelsberg telescope three decades ago. We compared 
these observations with the current spectra and found ve- 
locity shifts at some positions which correspond to an ac- 
celeration of the gas flow of V < 0.03 km s _1 yr _1 . 

The measured kinetic temperatures lie between 9 K and 
12 K for the starless and isolated star-forming sources, ex- 
cept for the fast rotator Do279P7 where Tki n ranges be- 
tween 12 K and 15 K. An increased value of Tki n in this 
case may be due to magnetic energy dissipation since mag- 
netic and nonthermal energy densities may be nearly equal 
in Do279P7: for n = 1.1 x 10 4 cm~ 3 , v = 0.4 km s _1 , and as- 



suming B ~ 20 ^G 3 we obtain B 2 /8ir w pv 2 /2 - 2 x 10~ n 
dyn cm~ 2 . However, to what extent the magnetic energy 
dissipation contributes to the gas heating requires further 
investigations. 

For the complex filamentary dark clouds, Tkin varies be- 
tween 9 K and 18 K and strongly fluctuates from point to 
point. We found no simple relationship between gas den- 
sity and kinetic temperature. High density condensations 
in the filamentary dark clouds have a large range of tem- 
peratures, most probably determined by radiation from em- 
bedded stars and dust. 

For the beam filling factor 77 = 1, we estimated the ex- 
citation temperature of an inversion doublet, the ammonia 
column densities, and H2 densities for those clouds where 
both the (1,1) and (2,2) transitions were detected. The typ- 
ical values of T ox lie between 3 K and 8 K, the total NH 3 
column densities between 0.3 x 10 15 cm~ 2 and 3 x 10 15 
cm~ 2 , and H2 volume densities between 0.3 x 10 4 cm~ 3 
and 4.2 x 10 4 cm" 3 . 

The study of Aquila dense cores will be continued, with 
the remaining ~ 100 targets scheduled for observations at 
the Effelsberg 100-m telescope in 2013. 
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Kawamura 01 (L694-2) 




150 100 50 -50 -100 -150 
R.A. Offset (arcscc) 



Fig. 1. (a): NH 3 (1,1) intensity map (fT MB dv) of KaOl. The 
starting point for the contour levels is 0.1 K km s _1 , the incre- 
ment is 0.7 K km s _1 between the first two and 0.8 K km s _1 
between the other contour levels. The crosses mark measured 
positions, (b): NH3(1,1) radial velocity (color map) structure. 
The velocity field is shown after subtraction of the mean radial 
velocity (Vlsr) = 9.576 km s _1 . The intensity peak is labeled in 
accord with Table 1. The Effelsberg beam size (HPBW) at 23.7 
GHz is 40". 



Dobashi 279 P6 




R.A. Offset (arcsec) 



Fig. 2. (a): Same as Fig. la but for the source Do279P6. The 
starting point for the contour levels is 0.5 K km s~ , the incre- 
ment is 0.5 K km s -1 between the first two and 1.0 K km s _1 
between the other contour levels, (b): NH3(1,1) radial velocity 
(color map) structure. The velocity field is shown after sub- 
traction of the mean radial velocity (V LS r) = 7.986 km s _1 . The 
intensity peaks are labeled in accord with Table 1. (c): NH"3(1,1) 
intensity map as in panel a and positions of YSO candidates ob- 
served with the Spitzer/IRAC in a star formation region Serpens 
cluster B (Harvey et al. 2006). Open circles of different col- 
ors mark Class I (red), intermediate type (yellow), and Class II 
(blue) YSOs. 
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Dobashi 279 P8 




R.A. Offset (arcscc) 



Fig. 3. (a): Same as Fig. la but for the source Do279P8. The 
contour levels start at 0.5 K km s _1 and increments are 0.5 
K km s _1 . (b): NH3(1,1) radial velocity (color map) structure. 
The velocity field is shown after subtraction of the mean radial 
velocity (Vlsr) = 8.553 km s _1 . The intensity peak is labeled in 
accord with Table 1. 



Dobashi 279 P12 




R.A. Offset (arcsec) 

Fig. 4. (a): Same as Fig. la but for the source Do279P12. The 
contour levels start at 0.5 K km s and increments are 1.5 
K km s _1 . (b): NH3(1,1) intensity map as in panel a and posi- 
tions of IR sources. Open circles mark the Herbig-Haro objects 
(blue), far- infrared continuum sources (red and yellow) (Davis 
et al. 1999), submm-continuum sources (magenta) detected by 
SCUBA (Di Francesco et al. 2008), and the triple radio source 
FIRS 1 (cyan) (Snell & Bally 1986). 
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Dobashi 279 P12 




300 250 200 150 100 50 -50 -100 
R.A. Offset (arcsec) 



Fig. 5. NH3(1,1) intensity (grey contours) and radial velocity ; ; 

(color map) structure in the source Do279P12. The velocity field -120 -bo -40 40 ao 120 ieo 

is shown after subtraction of the mean radial velocity (Vlsr ) = Ai (arcsec) 

8.099 km s . The contour levels are the same as in Fig. 4a. 
Panels A and B represent the velocity maps for the components 

A and B in the NH 3 profile, respectively (see Fig. B.4). The Fig. 6. Dots with 3ct error bars mark the differences between the 

white crosses mark offsets where the splitting of the NH 3 profiles radial velocities measured with the Effelsberg 100-m telescope 

is detected. The mean radial velocities of the components A in 1977-1981 (Ungerechts & Giisten 1984) and in the present 

and B are (Va) = 7.27 ± 0.11 km s _1 and (Vb) = 8.54 ± 0.08 dataset toward the same positions in Do279P12. The red points 

km s" 1 . In panel B, the intensity peaks are labeled in accord indicate deviations from zero which are larger than 3a. At two 

with Table 1. offsets, (Aq, A<5) = (80", 0") and (160", -120"), we detected 

double NH3 profiles (see Fig. B.4). 
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Dobashi 321 P2 




150 100 50 -50 -100-150-200 
R.A. Offset (arcsec) 



Fig. 7. (a): Same as Fig. la but for the source Do321P2. The 
contour levels start at 0.3 K km s _1 and increments are 1.0 
K km s . (b): NH3(1,1) radial velocity (color map) structure. 
The velocity field is shown after subtraction of the mean radial 
velocity (Vlsr) = 6.905 km s _1 . The intensity peaks are labeled 
in accord with Table 1. (c): NH3(1,1) intensity map as in panel a. 
Open circles of different colors mark Class I (red), intermediate 
type (yellow), Class II (blue), and Class III (cyan) YSOs taken 
from Tsitali et al. (2010). The SCUBA sources from Visser et 
al. (2002) are marked by magenta: SMM7 is a starless core, 
whereas SMM1 and SMM2 are protostars which coincide with 
IRAS 19180+1116 and IRAS 19180+1114, respectively. 



Kawamura 05 




40 -40 

R.A. Offset (arcsec) 



Fig. 8. (a): Same as Fig. la but for the source Ka05. The 
contour levels start at 0.4 K km s _1 and increments are 1.0 
K km s" 1 . The square indicates the position of the far-infrared 
source IRAS 19345+0727 coinciding with the 3.6 cm continuum 
source (Anglada et al. 1992). (b): NHa(l,l) radial velocity (color 
map) structure. The velocity field is shown after subtraction of 
the mean radial velocity (Vlsr) = 8.390 km s _1 . The intensity 
peak is labeled in accord with Table 1. 
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Dobashi 279 P7 




80 40 -40 -80 -120 -160 
R.A. Offset (arcscc) 



Fig. 9. (a): Same as Fig. la but for the source Do279P7. The 
contour levels start at 0.3 K km s _1 and increments are 0.6 
K km s _1 . (b): NH3(1,1) radial velocity (color map) structure. 
The velocity field is shown after subtraction of the mean radial 
velocity (Vlsr) = 6.715 km s _1 . The intensity peak is labeled in 
accord with Table 1. 



SS 3 




R.A. Offset (arcsec) 



Fig. 10. (a): Same as Fig. la but for the source SS3. The 
contour levels start at 0.0 K km s _1 and increments are 2.0 
K km s~ . (b): NH3(1,1) radial velocity (color map) structure. 
The velocity field is shown after subtraction of the mean radial 
velocity (Vlsr) = 7.542 km s . The white crosses mark offsets 
where the splitting of the NH3 profiles is detected (see Fig. B.8). 
The mean radial velocities of the components A and B are 
(Va) = 7.865 ± 0.005 km s" 1 and (V B ) = 6.856 ± 0.010 km s" 1 . 
The intensity peaks are labeled in accord with Table 1. (c): 
NH3 intensity map as in panel a and positions of the MAMBO 
millimeter sources taken from Table 1 in Maury et al. (2011). 
Red and yellow open circles mark starless and YSO candidates, 
respectively. 
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Table 1. Peak intensities of NH3(1,1) line emission toward Aquila targets 
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2.5(5) 
1.5(3) 
0.20(4) 
1.0(2) 

0.15(3 
0.10(2 

3.0(6 
2.7(6 
1.7(3 
1.6(3 
1.3(3 

2.5(5) 
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0.311(8) 



30-03-11 

26,29-03-11 
26,29-03-11 
26,29-03-11 
26,29-03-11 
26,29-03-11 

28-03-11 
28-03-11 
28-03-11 

13-01-12 
13-01-12 
12-01-12 

30-03-11 

12-01-12 

28-03-11 

28- 03-11 

29- 03-11 
29-03-11 

29-03-11 
29-03-11 
29-03-11 
29-03-11 
29-03-11 

12-01-12 

28-03-11 



J180447.8-043208 1 
Ser G3-G6 2 



BDN 31.57+5.37' 4 



BDN 28.65+3. 66 3 



BDN 46.22-1.34, 
L673, P77 3 



B335 4 , CB199 5 
B143 4 , L694-2 6 



Notes. "Greek letters mark peaks of ammonia emission indicated in Figs. 1-3, 5, 7-10. 

b The number in parentheses correspond to a la statistical error at the last digital position. c Linewidth (FWHM) . 
References. (1) Di Francesco et al. 2008; (2) Cohen & Kuhi 1979; (3) Dobashi et al. 2005; (4) Barnard 1927; 
(5) Clemens and Barvainis 1988; (6) Lee and Myers 1999. 
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Table 2. Molecular cores toward Aquila targets without detected NH3 emission in the (1,1) and (2,2) inversion transitions 
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Appendix A: The analysis of NH 3 (1,1) and (2,2) 
spectra 

In this Appendix we discuss the methods used for calculat- 
ing physical parameters from observations of the ammonia 
NH3 (J, K) = (1, 1) and (2,2) transitions. We give a consis- 
tent summary of the basic formulae relevant to the study 
of cold molecular cores. 

The recorded NH3 spectra were first converted into the 
T MB scale and a baseline was removed from each of them. 
The baseline was typically linear in the range covered by 
the ammonia hyperfine (hf ) structure lines but occasionally 
quadratic. After that the spectra of the same offset were 
co-added with weights inversely proportional to the mean 
variance of the noise per channel, rms -2 , estimated from 
the channels with no emission. 

The NH3 spectra were fitted to determine the total op- 
tical depth Ttot in the respective inversion transition, the 
LSR velocity of the line V^ SR , the intrinsic full width to 
half power linewidth Av for an individual hf component, 
and the amplitude, A (e.g., Ungerechts et al. 1980): 



T{v) = A 



-t(v) 



(A.l) 



The optical depth t(v) at a given radial velocity v is defined 

by 



t(v) 



Ttot ^V;exp{-2.77[(v- 

i=l 



V^ K )+ Vi } 2 /(Av) 2 }.(A.2) 



In this equation, the sum runs over the n hf components 
of the inversion transition (n = 18 and 21 for the NH3(1,1) 
and (2,2) inversion transitions, respectively), T{ is the the- 
oretical relative intensity of the zth hf line and Vi is its 
velocity separation from the fiducial frequency. The values 
of these parameters are given in, e.g., Kukolich (1967) and 
Rydbeck et al. (1977). 

In Eq.(A.2), the total optical depth, r to t, is the max- 
imum optical depth which an unsplit (1,1) or (2,2) line 
would have at the central frequency if the hf levels were 
populated at the same excitation temperature for both lines 
(1,1) and (2,2). Assuming that the line profile function has 
a Gaussian shape with the width Av (FWHM) and tak- 
ing into account that the statistical weights of the upper 
and lower levels of an inversion transitions are equal, one 
obtains 



Ttot (J, ^) = 



167TV / In 2 \ 1/2 K 2 
3hAv \ 7T J J(J + 
1 - cxp(-hv/k B T cx ) 
1 + ex.p(-hv/k B T eK ) ' 



-N(J,K) 



(A.3) 



where h and k B are Planck's and Boltzmann's constants, 
respectively, fj, is the dipole moment, = 1.4719 Debye 
(JPL catalog 4 ), N(J,K) is the inversion state column den- 
sity — the sum of the column densities of the upper and 
lower levels of an inversion doublet, v is the rest frequency 
of the inversion transition, and T cx is the excitation tem- 
perature which characterizes the NH 3 population across the 
(J, K ) inversion doublet. 



http://spec.jpl.nasa.gov/ 



In Eq.(A.l), the amplitude A can be expressed in terms 
of the beam filling factor, r\ = r2 cloud /O boam (solid an- 
gle il = 7r# 2 /(41n2), where 9 is the angular diameter at 
FWHP), and the excitation temperature, T cx , as 



A = rj[J v {T m )- MT hg )] , 



(A.4) 



where Tbg = 2.7 K is the black body background radiation 
temperature, and J V (T) is defined by 



J V (T) = — e ^/ feBT -l 



r 



(A.5) 



For optically thin lines (e.g., m -C 1), Eq.(A.l) degen- 
erates into 

T(v) = A-t{v) 

n 

= {Am) ■ J2 n cxp[-2.77(A\yAi;) 2 ] , (A.6) 

i=l 

and it is not possible to determine A and m separately. 
Here AVi = v- V LSR + v t . 

For some molecular cores we observe two-component 

NH3 profiles. The line parameters {Aj, t$ , , V^siL Avj} 
for both j = 1 and j = 2 components were determined in 
this case by fitting the function 



T(v)=T 1 {v)+T 2 (v). 



(A.7) 



This linear form assumes that two clouds in the beam have 
the same filling factor r\ and the same T ex for both lines 
(1,1) and (2,2). 

At low spectral resolutions, when the channel spacing, 
A c h, is comparable to the linewidth Av, the derived line 
parameters are affected by the channel profile </>(i>)- In this 
case, the recorded spectrum, T'(v), was considered as a 
convolution of the true spectrum T(v) and 4>{v) 



7» 



/t«,) 



4>{v — v ) dv . 



(A.8) 



We approximate <j>(v) by a Gaussian function with the 
width FWHM defined by a particular backend setting. 

The calculated synthetic spectrum which is defined as 
T syn (v) = T(v) if A^ « Av, or T syn (v) = T'(v) if A 



^ch 



Av, was fitted to the observed spectrum T b s (v) through 
minimization of the x 2 function 



X 



[T syn (i;)-r obs (i;)] 2 /rnis 2 



(A.9) 



in the 5-dimcnsional space of the following model param- 
eters: A, Tn,T22, V LSR , and Aw. Attempting to exclude at 
the 3<T level all noise-only channels in between the hf lines, 
the sum runs over the channels with observed NH3 emission 
from all hf components. The minimum of \ 2 was computed 
using the simplex method (e.g., Press et al. 1992). 

Equation (A.9) was also used to estimate the formal 
errors of the model parameters through calculation of the 
covariance matrix at the minimum of x 2 . Since the un- 
certainty in the amplitude scale calibration was ^15-20% 
(Sect. 2), and there were no noticeable correlations between 
the sequential channels in our datasets, we did not correct 
the calculated errors by an additional factor as described, 
e.g., in Rosolowsky et al. (2008). 



18 



S. A. Levshakov et al.: NH3 observations in the Aquila rift cloud complex 



Given the estimate of the amplitude A, the excitation 
temperature T ox can be calculated from Eq.(A.4): 



To 



ln[l + T QV /(A + r,J hg )} 



(A.10) 



where T = hv/kv, and J bg = Ju{T hg ). At T hg = 2.725 K 
and NH 3 -inversion frequencies, we have Jb g = 2.20 K and 
T a = 1.14 K. 

In this equation the only unknown parameter is the 
filling factor, 77. If the source is completely resolved, then 
77 = 1. Otherwise, Eq. (A. 10) gives a lower bound on T cx 
corresponding to rj = 1, whereas decreasing r\ drives T cx to- 
wards its upper bound at T cx = Tkin, which holds for LTE. 
In case of emission, the radiation temperature within the 
line is larger than the 2.7 K background and the excitation 
may be either radiative or collisional. At gas densities larger 
than the critical density 5 and for optically thin lines the lo- 
cal density of line photons is negligible as compared to the 
background radiation field, and the lower metastable states 
in NH 3 are mainly populated via collisions with molecular 
hydrogen 6 , i.e., the stimulated emission can be neglected. 
For low collisional excitation, n gas < n cr , the line radia- 
tion controls T ex for both optically thin and thick regimes 
(Kegel 1976). NH 3 (1,1) emission usually arises from regions 
with n gas > 10 4 cm~ 3 (e.g., Ho & Townes 1983). It means 
that the regime of low collisional excitation does not domi- 
nate. We also do not consider anomalous excitation (maser 
emission), leading to T cx >• Tb g , since masing is not ob- 
served in our data sample. Thus, the use of T cx — 21cin as a 
formal upper limit for the excitation temperature is legiti- 
mate. By this we can restrict the unknown filling factor for 
an unresolved source in the interval 



?7min ^ T] — ?7max ; 

and T cx within the boundaries 



lr,=l 



rrmax I 



(A.11) 



(A.12) 



From a separate analysis of the (1,1) and (2,2) lines 
with low noise level, we found that both amplitudes An and 
A22 are equal within the observational errors, which implies 
that T cx (l, 1) s=s T cx (2, 2). Albeit radiative transfer calcula- 
tions show that the excitation temperatures T ox (l, 1) and 
T ex (2,2) may differ by about 20% (Stutzki & Winnewisser 
1985), the assumption of their equivalence seems to be suf- 
ficiently good in our case. 

The second temperature which describes the NH3 pop- 
ulation is the rotational temperature, T rot , characterizing 
the population of energy levels with different (J, if). As 
mentioned above, the population of the lower metastable 
inversion doublets is determined by collisions with H2 and, 
thus, is regulated by the kinetic temperature, Tkin- The pop- 
ulation ratio between the (1,1) and (2,2) states is defined 
as 



12 



^22 _ ff22 c „ / A£; 
Nu 511 \ k B T I0 t 



-•exp 



41. 5K 



trot 



. (A.13) 



5 For NH 3 (1,1), n cr = 3.90 x 10 3 cm" 3 at T kin = 10 K (Maret 
et al. 2009). 

6 The abundance ratio [Hhj/pHe] ~ 5 and the NH3-H2 collision 
rates are a factor of 3 larger than the NH3-He rates (Machin & 
Roueff 2005). 



Assuming that both transitions trace the same volume 
of gas and their linewidths Av are equal, one finds from 
Eqs.(A.3) and (A.13) that 



T rot = -41.5/ In 



20 



T22 

Til 



(A.14) 



We note that since radiative transitions between the differ- 
ent if-laddcrs of the (J, K) levels are forbidden and that 
the unknown filling factor 77 is canceled in the line intensity 
ratio of the NH3 (1,1) and (2,2) transitions, the estimate of 
the rotational temperature is more reliable than the exci- 
tation temperature. 

For a two-level system with the kinetic temperature less 
than the energy gap between the (1,1) and (2,2) states 
(A_Ei2 = 41.5 K), the rotational temperature can be re- 
lated to the kinetic temperature through detailed balance 
arguments (e.g., Ragan et al. 2011): 



Tkin — T rot /[l — 



AE 



12 



ln(l + 1.1 c- 15 - 7 / 71 ™ 1 )] . (A.15) 



The beam-averaged column density N JK (in cm -2 ) can 
be calculated from Eqs.(A.3) and (A. 4) using the estimated 
values of the amplitude A, the total optical depth t jk , and 
the linewidth Aw: 



N JK = (.IK ■ Av ■ Tjf 



A + T]( Jbg + 0.5T ) 



(A.16) 



where Cn = 1-3850 x 10 13 , ( 22 = 1.0375 x 10 13 , and Ai; 
is in km s _1 (for detail, see Ungerechts et al. 1980). As 
mentioned above, An ~ A22, and for both transitions T = 
1.14 K, and Jbg = 2.20 K. Substituting these numerical 
values in Eq.(A.16), one obtains 



N n = 3.84 x 10 13 • Av ■ r n • (1 + -4/2.77?7) • 



(A.17) 



This gives us lower and upper boundaries of Nu if 77 e 

[?7min, 1] : 



1,7=1 



< Nn < N? 



<V = Vn 



(A.18) 



The uncertainty interval AN n = N^ x - A^f n trans- 
forms directly into the estimate of the total NH3 column 
density: 



AT(NH 3 ) = N n ■ [(1/3) exp(23.2/T rot ) + 1 + 
+(5/3)exp(-41.5/T rot ) + 
+ (14/3)exp(-105.2/T rot ) + ...] . 



(A.19) 



which assumes that the relative population of all metastable 
levels of both ortho-NH3 (K = 3), which is not observable, 
and para-NH3 (K — 1,2) is governed by the rotational tem- 
perature of the system at thermal equilibrium (Winnewisser 
et al. 1979). Substituting iV™ in and iVft 8 * from (A.18) into 
(A.19), one finds: 



iV(NH 3 ) r 



=1 < A^NHg) < N(NH 3 ) 



(A.20) 



The detailed balance calculations also provide a relation 
between the gas density, gas kinetic temperature, and the 
excitation temperature (e.g., Ho & Townes 1983): 



»(H 2 ) = ^ 



Jv{T cyL ) — Jy(Tbg) 



Jf (^kin) J v (^cx) 



1 + 



J/ (Tkin) 

hv/k B 



(A.21) 
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where A is the Einstein A-coefficicnt, and C is the rate 
coefficient for collisional de-excitation. For a typical kinetic 
temperature in the dense molecular cores of <~ f 0— 20 K, the 
collision coefficient is ~ 4 x f0 _7 s _1 [n(H 2 )/T0 4 cm -3 ]. The 
value of A for the inversion transition (f,f) is f.67 x 10~ 7 
s _1 . However, the gas density calculated using Eq.(A.21) 
may be significantly underestimated if the beam is not filled 
uniformly. Moreover, in case of T cx = Tki n , Eq.(A.2f) is 
invalid and, hence, n(H 2 ) is to be calculated by different 
methods (see, e.g., Hildcbrand f983; Pandian et al. 2012). 
We use (A. 21) just to set a lower bound on the gas density, 
n(H 2 ) min , at 77 = 1. 

Given a fractional NH 3 abundance, an upper bound 
on the gas density may be estimated from the deduced 
Af(NH 3 ) max assuming that the ammonia emission traces 
the real distribution of the gas density 7 . If the source is 
unresolved, its diameter, d, and the beam filling factor are 
related as 



d = 9 m D^j 



(A.22) 



where 6* m is the beam angular diameters (FWHP), and D 
the distance of the source. The highest gas density is ob- 
tained at the smallest diameter, d mlD = 9 ul D y /i r i m i n : 



n(H 2 ) r 



iV(NH 3 ) r 



X-d n 



(A.23) 



To conclude, we note that in some cases, when Tkin 
varies slowly within an ammonia clump, the NH 3 column 
density can be deduced from the (1,1) line alone assuming 
the same average Tkin in the core as in the outskirts of a 
cloud. This allows us to extend the radial gas distribution 
to zones with unobservable (2,2) emission (Morgan et al. 
2013). 



where X = [NH 3 ]/[H 2 ] is a given abundance ratio. In this 
equation the unknown distance D can be found from the 
requirement that both values of the lowest gas densities 
calculated at 77 = 1 from Eq.(A.21) and from -/V(NH 3 ) m j n 
at rf max = m D be equal. This gives 



n(H 2 ) r 



n(H 2 ) jV(NH 3 ) max 

min 



(A.24) 



For the unresolved source its axis — major, 9 a , or mi- 
nor, 9b — is less than 9 m , rj < 1, and only limiting values 
can be obtained for T ox , Njk, N(NB.s), and n(H 2 ). If the 
source is resolved, i.e., #b > 9 m , then rj — 1 and these 
physical parameters are directly defined. In the latter case 
we calculate the deconvolved values of 9 a and 9^ and their 
geometrical mean 



= V^a • 9 h 



(A.25) 



which is used as a formal estimate of the source angular 
diameter. 

A similar procedure is applied to calculate the virial 
mass of the ammonia core through the linewidth, Aw. If 
the observed value Aw is of order of the spectral resolution 
A sp (FWHM), then only an upper limit on M v ; r is defined. 
Otherwise, if Aw > A sp , then the deconvolved value of the 
linewidth, Av, and the core radius r give the virial mass 
(e.g., Lemme et al. 1996) 



250 • Av 2 ■ r 



(A.26) 



where Av is in km s \ r in pc, and M v i r in solar masses 
M . 



Ammonia and submillimeter maps of dense cores in 
Ophiuchus show a close correlation between the large-scale dis- 
tributions of the NHa(l,l) integrated intensity and the 850 /im 
continuum emission (Friesen et al. 2009). 
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Appendix B: Ammonia spectra toward the Aquila 
rift cloud complex and derived physical 
parameters 

The observed spectra of the NH 3 (1,1) and (2,2) transitions 
detected at the peak positions of ammonia emission toward 
each core are shown in Figs. B.1-B.9. The measured phys- 
ical parameters are listed in Tables B.1-B.4. The distribu- 
tions of the physical parameters for three most abundant 
cores are presented in Figs. B.10-B.12. They demonstrate 
variations of the kinetic temperature, Tkin, the excitation 
temperature, T G x, the linewidth, Av (FWHM), the total 
optical depth, rn, the gas density, hh 2 , and the total am- 
monia column density, iV(H2), across the mapped area of 
the core. 




Fig. B.l. Ammonia (1,1) and (2,2) spectra (blue) toward the 
source KaOl. The channel spacing is 0.015 km s _1 , the spectral 
resolution FWHM = 0.024 km s" 1 . The red curves show the fit 
of a single-component Gaussian model to the original data. The 
residuals between the observed and model spectra are shown by 
black color. 




NH 3 (2,S) 



- Do279P8 NH.,0,1) offset (0,0) 



- NH,(Z,2) 




-10 -5 5 10 15 20 25 5 10 

V LSB (km S">) V LSH (km s"«) 

Fig. B.3. Same as Fig. B.l but for the source Do279P8 and the 
channel spacing 0.077 km s" 1 (FWHM = 0.123 km s _1 ). The 
red curves show the fit of a single-component Gaussian model 
to the NHs(l,l) original data and the upper limit on NHg(2,2). 



Bo^7BPTl2 NH 3 (1 ,l')A ' offset ' (20o! -1 20) 




yv u^^v 




Fig. B.4. Same as Fig. B.l but for the source Do279P12. The 
ammonia spectra have double structure. The upper panels show 
the high-resolution spectra (channel spacing 0.015 km s _1 ), and 
the other panels — low-resolution spectra (channel spacing 0.077 
km s _1 ). Two NH3 components are marked by ticks and labeled 
by letters A and B in the upper panels. 




Fig. B.5. Same as Fig. B.l but for the source Do321P2. 



-10 -5 5 10 15 20 25 5 10 

V lsb ( km s ~") V lse ( km s ~') 



Fig. B.2. Same as Fig. B.l but for the source Do279P6. 
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Fig. B.6. Same as Fig. B.l but for the source Ka05 and the 
channel spacing 0.038 km s^ 1 (FWHM = 0.044 km s _1 ). 



z 


-2 

-, 2 
6 

_S 

-2 
2 



Do279P? NII 3 ( 1 , 1 ) 



offset (-40,40) 



1 1 1 1 1 1 1 1 1 1 1 1 1 

offset ( 80, Q) 



I I I I I I I I I I I I 
offset ( 80,40) :: 



offset. (-80.80) 



NH 3 (2,2) 



-15 -10 -5 5 10 15 20 25 5 10 

V LSR ( km S_1 ) V LSR ( klT1 S_1 ) 



Fig. B.7. Same as Fig. B.6 but for the source Do279P7. 



SS3 NH 3 (1,1) 



8 3 




offset (8Q, 160) 



— H — H — 



offset (0,160) 



I I 1 1 1 1 1 1 1 1 1 



offset (0,80) 




NH 3 (2,2) 



Fig. B.8. Same as Fig. B.l but for the source SS3. The chan- 
nel spacing is 0.77 km s _1 and 0.038 km s _1 at two upper and 
three lower panels, respectively. The corresponding spectral res- 
olutions are 0.895 km s _1 and 0.044 km s" 1 (FWHM). The red 
curves show the fit of a single-component (two upper panels) 
and a double-component (three lower panels) Gaussian model 
to the original NH3 data. Two components of the NH3 emission 
are marked by ticks and labeled by letters A and B. 



-Do243 P2 NH 3 (1,1) offset. (4 0,0) 



-Do279 P18 NH 3 (l,l) offset (0,0) 



f-H — M ill 



II 



Do279 P13 NH 3 (1,1) 1 offset (-40,0) 
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Do321 PI NH 3 (1,1) offset (-40,0) 




-Do32i PI NH 3 (1,1) 



offset (0,40) - 



ft 



-10 -5 5 10 15 20 25 
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Fig. B.9. Same as Fig. B.3 but for the sources Do243P2, 
Do279P18, Do279P13, and Do321Pl. 
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Fig. B.10. The spatial distributions of the physical parameters measured in Do279P6 from the ammonia inversion lines NH3(1,1) 
and (2,2). The corresponding numerical values are listed in Table B.2. 
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Dobashi 279 P12 




Fig. B.ll. The spatial distributions of the physical parameters measured in Do279P12 from the ammonia inversion lines NH3(1,1) 
and (2,2). The corresponding numerical values are listed in Table B.3. 
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Fig. B.12. The spatial distributions of the physical parameters measured in SS3 from the ammonia inversion lines NH3(1,1) and 
(2,2). The corresponding numerical values are listed in Table B.4. 
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Table B.l. Observed parameters of the NH3(1,1) and (2,2) lines and calculated model parameters for Kawamura 01 and 05, 
Dobashi 279 P7 and P8, and Dobashi 321 P2 



Peak 


Aa AS 


Vlsr 


Av 


T 


Tex 


T 


-^kin 


riH 2 /10 4 




T 22 


iVtot/10 15 


Id. 


(") (") 


(km s- 1 ) 


(km s- 1 ) 


(K) 


(K) 


(K) 


(K) 


(cm" 3 ) 






(cm^ 2 ) 










rvaUl 


















0,0 


9.60 


0.32 


1.5 


4.3 


8.7 


9.1 


1.2 


10.9 


0.2 


1.2 


a 


0,-40 


9.55 


0.32 


2.5 


5.5 


9.2 


9.6 


2.5 


10.5 


0.3 


1.3 




+40^-40 


9.60 


0.30 


1.4 


4.8 


< 8.5 






4.1 


< 0.1 


> 0.5 




+40,-80 


9.54 


0.30 


0.8 


4.1 


< 9.6 






3.2 


< 0.1 


> 0.3 




-40,0 


9.62 


0.23 


1.7 


4.7 


< 8.9 






8.2 


< 0.2 


> 1.0 




-40,-40 


9.60 


0.27 


1.3 


4.3 


< 9.7 






6.6 


< 0.2 


> 0.5 










Ka05 
















a 


0,0 


8.40 


0.37 


2.5 


5.6 


10.2 


10.7 


2.3 


5.6 


0.3 


0.8 




0,-40 


8.40 


0.33 


2.1 


5.1 


10.8 


11.5 


1.6 


8.0 


0.4 


0.7 










Do279P7 














a 


-80,40 


6.04 


0.75 


1.5 


6.4 


12.6 


13.7 


1.1 


4.3 


0.4 


0.6 




-40,40 


6.14 


0.72 


1.4 


4.5 


13.3 


14.7 


1.0 


3.8 


0.4 


0.5 




-80,0 


6.59 


0.55 


1.4 


5.1 


13.5 


14.9 


1.4 


2.4 


0.3 


0.3 




-80,80 


5.91 


0.73 


0.5 


3.3 


11.5 


12.4 


0.3 


4.4 


0.3 


0.5 










Do279P8 














a 


0,0 


8.52 


0.27 


2.5 


5.7 


< 8.0 






8.0 


< 0.1 


> 1.1 




0,-40 


8.57 


0.30 


2.4 


5.7 


< 7.7 






6.9 


< 0.1 


> 1-2 










Do321P2 














a 


-80,120 


6.88 


0.48 


3.0 


6.8 


11.9 


12.9 


3.3 


4.0 


0.3 


0.6 


P 


-80,40 


7.05 


0.43 


2.7 


6.7 


11.4 


12.2 


3.4 


3.5 


0.2 


0.5 


7 


-120,-40 


6.93 


0.38 


1.7 


5.5 


< 9.1 






3.0 


< 0.1 


> 0.5 


S 


-80,240 


7.12 


0.48 


1.6 


4.7 


< 9.7 






4.8 


< 0.1 


> 0.7 


e 


80,-40 


6.66 


0.73 


1.3 


< 5.7 















Notes. An is the linewidth (FWHM); m, r 22 are defined in Eqs. (A.2) and (A. 3), and JV to t — in Eq. (A. 19). 



Table B.2. Observed parameters of the NHs(l,l) and (2,2) lines and calculated model parameters for Do279P6 



Peak 


Aa AS 


Vlsr 


Av 




T 

- 1 ex 


Trot 


Tkin 


n H2 /10 4 


Til 


T~22 


iVtot/10 15 


Id. 


(") (") 


(km s" 1 ) 


(km s" 1 ) 


(K) 


(K) 


(K) 


(K) 


(cm" 3 ) 






(cm- 2 ) 


a 


0,0 


7.75 


0.75 


2.9 


5.7 


12.8 


14.1 


1.9 


5.7 


0.5 


1.1 


P 


40,40 


8.04 


0.63 


2.1 


5.3 


10.8 


11.5 


1.8 


6.0 


0.3 


1.1 


7 


-160,-80 


7.69 


0.96 


1.8 


5.8 


12.9 


14.2 


2.0 


2.2 


0.2 


0.5 


8 


-120,-40 


7.89 


0.68 


1.7 


5.6 


10.6 


11.3 


2.2 


2.7 


0.1 


0.6 


e 


80,80 


7.81 


0.43 


1.3 


4.1 


9.5 


9.9 


0.8 


10.9 


0.3 


1.3 




40,0 


7.82 


0.90 


1.2 


4.5 


13.6 


15.1 


1.0 


2.7 


0.3 


0.5 




0,40 


7.97 


1.08 


1.5 


5.2 


14.2 


16.0 


1.4 


2.4 


0.3 


0.5 




0,-40 


7.81 


0.73 


1.6 


4.7 


14.2 


15.5 


1.1 


4.4 


0.5 


0.6 




-40,-40 


7.95 


0.77 


1.4 


5.4 


15.3 


17.5 


1.5 


2.0 


0.3 


0.3 




-40,0 


7.81 


0.95 


1.2 


4.2 


11.6 


12.5 


0.9 


4.5 


0.3 


0.9 




-80,0 


8.08 


0.38 


1.2 


5.2 


13.2 


14.5 


1.4 


2.4 


0.2 


0.2 




80,40 


7.78 


0.40 


0.9 


3.9 


12.2 


13.2 


0.6 


6.0 


0.4 


0.4 




-80,-40 


8.15 


0.73 


1.1 


5.1 


12.8 


14.0 


1.4 


1.6 


0.1 


0.3 




-160,-40 


7.78 


1.00 


1.5 


5.1 


12.2 


13.2 


1.5 


2.3 


0.2 


0.5 




-200,-80 


7.96 


0.83 


2.2 


5.5 


11.9 


12.9 


1.9 


4.2 


0.3 


0.9 




-200,-40 


7.98 


0.85 


1.2 


4.5 


10.1 


10.6 


1.1 


3.3 


0.1 


0.8 




-200,-120 


8.19 


0.68 


2.2 


5.3 


11.3 


12.1 


1.7 


6.7 


0.4 


1.2 


Notes. Aw is the linewidth (FWHM); m, 


T22 are 


defined 


in Eqs. (A.2) and (A.3), and N tot - 


in Eq. (A. 19). 
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Table B.3. Observed parameters of the NH3(1,1) and (2,2) lines and calculated model parameters for Do279P12 



Peak 


Aa AS 




Aw 






Trot 


Tkin 


n H2 /10 4 
(cm^ 3 ) 


Til 


T~22 


iVtot/10 15 


Id. 


(") (") 


(km s" 1 ) 


(km s" 1 ) 


(K) 


(K) 


(K) 


(K) 






(cm" 2 ) 


a 


160,-40 


7.59 


0.90 


3.2 


7.6 


15.6 


17.8 


3.2 


2.8 


0.4 


0.7 





200,-120° 


A 6.94 


0.75 


3.1 


7.4 


13.5 


14.9 


3.5 


3.9 


0.4 


0.8 




B 8.21 


0.73 


1.2 


5.3 


14.1 


15.8 


1.5 


1.7 


0.2 


0.3 




200,-200 a 


A 7.26 


0.62 


2.3 


7.5 


11.7 


12.6 


4.5 


1.9 


0.1 


0.4 




B 8.58 


0.77 


2.0 


5.8 


12.6 


13.7 


2.1 


2.8 


0.2 


0.6 


7 


40,40 


8.29 


0.90 


2.9 


6.4 


12.2 


13.2 


2.8 


4.1 


0.3 


1.1 


0,0 


8.23 


1.22 


1.3 


4.9 


14.6 


16.5 


1.1 


2.3 


0.3 


0.5 




40,0 


8.42 


1.15 


1.7 


5.1 


12.8 


14.0 


1.4 


3.0 


0.3 


0.8 




80,0° 


A 7.89 


0.68 


1.3 


6.9 


16.4 


19.0 


2.5 


1.0 


0.2 


0.2 




B 8.92 


0.83 


0.6 


5.1 


14.1 


15.8 


1.3 


0.7 


0.1 


0.1 




0,80 


8.44 


0.98 


1.8 


5.3 


15.8 


18.2 


1.4 


3.1 


0.5 


0.6 




40,80 


8.35 


0.77 


2.6 


6.0 


12.9 


14.1 


2.2 


4.5 


0.4 


0.9 




40,120 


8.42 


0.68 


2.8 


6.5 


12.8 


14.0 


2.6 


4.0 


0.4 


0.8 




40,-40 


7.91 


0.77 


1.3 


4.9 


13.8 


15.4 


1.2 


2.4 


0.3 


0.4 




0,120 


8.56 


1.03 


2.0 


5.6 


14.7 


16.6 


1.7 


3.0 


0.4 


0.7 




80,80 


8.19 


0.97 


1.3 


4.3 


13.9 


15.4 


0.8 


4.3 


0.5 


0.7 




80,40 


8.51 


0.80 


2.2 


5.5 


11.8 


12.8 


1.9 


4.5 


0.3 


0.9 




120,-40 


8.16 


0.73 


1.6 


5.0 


13.8 


15.4 


1.2 


3.6 


0.4 


0.5 




120,-80 


8.15 


1.17 


1.8 


7.0 


12.3 


13.4 


3.4 


1.3 


0.1 


0.5 




160,-80 


7.67 


1.25 


2.4 


6.4 


14.0 


15.6 


2.4 


2.6 


0.3 


0.8 




160,-120° 


A 7.02 


0.68 


1.6 


7.0 


14.5 


16.2 


2.9 


1.3 


0.2 


0.2 






B 8.39 


0.90 


1.5 


5.6 


14.8 


16.8 


1.6 


1.9 


0.3 


0.4 




160,-160° 


A 7.40 


1.5 


1.0 










0.3 








B 8.64 


0.7 


1.3 










1.3 








160,-200 


A 7.18 


0.3 


1.1 










1.0 








B 8.59 


1.0 


1.0 










1.4 








240,40 


7.80 


0.73 


1.6 


5^6 


13.6 


15.1 


1.7 


2.3 


6.2 


0^4 




200,-40 


7.88 


0.87 


2.1 


10.2 


15.2 


17.3 


7.0 


0.8 


0.1 


0.3 




200,-80 


7.06 


1.28 


1.9 


6.5 


14.1 


15.7 


2.5 


1.5 


0.2 


0.5 




200,-160° 


A 7.26 


0.62 


2.3 


7.5 


11.7 


12.6 


4.5 


1.8 


0.1 


0.4 




B 8.58 


0.77 


2.0 


5.8 


12.6 


13.7 


2.1 


2.8 


0.2 


0.6 




200,-240 


8.29 


0.83 


1.8 


5.2 


12.8 


14.1 


1.5 


3.4 


0.3 


0.6 




200,-280 


8.42 


0.60 


1.9 


5.6 


11.3 


12.1 


1.9 


3.5 


0.2 


0.6 



Notes. "Partly resolved velocity components marked by letters A and B are shown in Fig. B.4. 

Av is the linewidth (FWHM); tii, r 22 are defined in Eqs. (A.2) and (A. 3), and JVtot — in Eq. (A. 19). 
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Table B.4. Observed parameters of the NH3(1,1) and (2,2) lines and calculated model parameters for SS3 
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Notes. "Partly resolved velocity components marked by letters A and B are shown in Fig. B.8. 

Av is the linewidth (FWHM); m, r 22 are defined in Eqs. (A.2) and (A.3), and Atot — in Eq. (A. 19). 



